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We give a very brief overview of collective effects in neutrino oscillations in core collapse supernovae where 
refractive effects of neutrinos on themselves can considerably modify flavor oscillations, with possible repercussions 
for future supernova neutrino detection. We discuss synchronized and bipolar oscillations, the role of energy and 
angular neutrino modes, as well as three-flavor effects. We close with a short summary and some open questions. 



1. Introduction 



sity matrices, 



Neutrinos interact via the weak interactions 
and can thus cause a refractive effect on each 
other. Under most circumstances, the resulting 
self-interaction potential is much smaller than the 
vacuum oscillation term or the potential induced 
by ordinary matter. However, in the early cool- 
ing phase of core collapse supernovae, the den- 
sity of neutrinos streaming off the hot nascent 
neutron star is sufhcicntly high to cause non- 
linear phenomena that can have practical im- 
portance. This was realized only recently in a 
series of papers of which we here cite only a 
few |ll:jl3l4l5l6l7l8l9ll0lllll2ll3lli] . 

To describe collective effects it is useful to de- 
scribe the neutrinos with flavor density matrices 
for each momentum mode. They are defined in 
terms of the annihilation operators a; for neutri- 
nos and (ij for antineutrinos in a given momentum 
mode p |15I16I17| . 
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Ip and {Q-p)ij = (a] ai)-p . (1 



their equations of motion can be written as a com- 
mutator of an effective Hamiltonian with the den- 



dtQp 



n„ + V2GfL- 
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Here, Gp is the Fermi constant, in the mass 
basis the matrix of vacuum oscillation fre- 
quencies for relativistic neutrinos is f2p = 
diag(m^, TO2, m§)/2p with p — |p|, and in the fla- 
vor basis L = dia,g(ne,n^,nT-), where n(,,n^,nT 
are the charged lepton number densities (parti- 
cle minus antiparticle density) responsible for the 
matter induced potential. When neutrino self- 
interactions are negligible, a compensation of the 
vacuum and matter induced terms in the diag- 
onal part of the Hamiltonian in Eq. ([2|) leads 
to the well known Mikheyev-Smirnov-Wolfenstein 
(MSW) effect jl8ll9| . An equation analogous to 
Eq. ([2]) holds for antineutrinos with the sign of 
the vacuum term Qp reversed. The last term 
in Eq. ^ describes forward scattering on the 
neutrino background where the isotropic den- 
sity term can lead to self-maintained coherence, 
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whereas the flux term proportional to the co- 
sine of the angle 6pq between momenta p and 
q can lead to self-induced decoherence. The flux 
of charged leptons has been assumed to vanish. 

In the two-flavor case it is convenient to 
parametrize vacuum, matter and self-interaction 
terms by the frequencies 



A(r) 



2p ' 
V2GFne(r); 

V2Gf [i^p, (r) - Fp^. (r)] (1 - cos 9^ 



(3) 



with Fi are the fluxes of the relevant neutrino 
species i, and stands for any one of z^^, 

or Vt. The oscillation phenomena discussed in 
the following depend crucially on the relative size 
of these three frequencies. 

2. Synchronized and Bipolar Oscillations 

Since muon and tau neutrinos behave equally 
in core collapse supernovae up to an effective /ir 
potential briefly discussed in Sect. [51 neutrino os- 
cillations can be described in two-flavor approx- 
imation. The oscillations are then driven by the 
atmospheric mass squared difference, |Am^jjj^| = 
|m§ - ml\ ^ 2.40+°;}^ x lO^^eV^ and the an- 



gle sin 6*13 < 0.04 



For these parameters. 



0.4 km^^ for a typical neutrino energy of 
15 MeV. Both normal and inverted hierarchy are 
still allowed for this case. 

Considering only a single two-flavor mode p 
one can expand into Pauli matrices cr , Q = {f + 
P • (t)/2, analogously for g, L — {uq + rigL • <t)/2, 
and u = [u;o + |Ato2|B • cr/(2p)] /2. With the 
vacuum mixing angle 0, B = (sin 26*, 0, — cos 26*), 
the radial flavor evolution (we use units in which 
the speed of light is unity) can be described by [4] 

drP = [+ujB + A(r)L + fi{r) (P - P)] x P , (4) 

and an analogous equation for P with the sign of 
Lu reversed. It is useful to define the vectors 
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varying /i, the equations of motion (now written 
in terms of time derivatives) are 



Q = ^DxQ; D = cjBxQ, 
which implies 



IQI — const ; 

D 



(7 = D x Q = const 

q X q 



o-q. 



(6) 



(7) 



These are the equations for a spinning top of spin 
cr, angular momentum (q x q)//i and moment of 
inertia / = mP — fJ.^^. Its energy then has a 
kinetic and a potential part, 

E = Fkin + Spot - |d2+C^(B-Q + Q) 

= ^ + ^a^+u;{B-Q + Q). (8) 

This description is a good approximation in the 
regime where neutrinos are freely streaming. The 
supernova neutrino (number) fluxes are thought 
to obey the hierarchy Fi,^ > Fj^^ > F^^. One 
can then adopt initial conditions at the neutrino 
sphere are typically P = (0, 0, 1+e), P = (0, 0, 1), 
where the asymmetry parameter in terms of the 
fluxes of different neutrino flavors is 



Fr. 



F„. 



1 



F,, - Fp^ 

Fr,. - Fr,^ 



(9) 



In the absence of matter, A = 0, and for slowly 



Eq. ^ then implies that the "flavor pendulum" 
is in a stable initial position for the normal hier- 
archy, 6^1, whereas for the inverted hierarchy 
9 = t:/2 — 9<^1 the pendulum is initially in a 
maximum energy state. We will focus on this case 
in the following. 

As long as > 2w/(l — \/l + e)^ the kinetic 
(self-interaction) term in Eq. ([8]) dominates and 
the oscillations are synchronized with a common 
frequency Wsynch = (2 -I- e)uj/£ around B. For 
u> < fj. < 2uj/{l — vT+T)^, bipolar oscillations 
take place with a frequency k = {2ujfi/ (1 + e))^/^ 
in which Q swings between its initial position and 
a position in which it is parallel to — B. This cor- 
responds to a collective transition of Vei^e to Vx^x 
pairs, with a rate greatly speeded up compared 
to ordinary pair annihilation [T]. Furthermore, 
in the absence of matter D • B = (P — P) • B 
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is strictly conserved ^ , whereas in dense matter 
(P — P) • L is approximately conserved. For small 
effective mixing angle, one thus has ^ Pz + 
Finally, when /i < w at large radii, vacuum oscil- 
lations ensue or, in the presence of matter with 
X ^ u! the usual MSW effects can occur. 

In the limit of small vacuum mixing angle, 9 = 
7r/2 — <C 1, and for constant fj, and A, the time 
scale for bipolar conversion is 



'^bipola 



-K^i In 



(At2 + A2)l/2 



(10) 



The main effect of the matter term is thus to 
decrease the effective mixing angle to 6k/{k^ + 
A^)-'^/^ gj^(\ i\ms delay the onset of bipolar oscil- 
lations. For smaller mixing angle it thus takes 
longer to tip over, the bipolar transition is less 
adiabatic and the nutation amplitude is larger. 
Since matter decreases the effective mixing angle, 
a high matter density also leads to a later onset 
of the bipolar transition and to a larger nutation 
amplitude 6J. 

For varying neutrino density and in the limit of 
small vacuum mixing angle, angular momentum 
and adiabatic energy conservation in the bipo- 
lar regime gives for the electron neutrino survival 
probability P{ve ~^ Ve) = \{l + Pz) oc 

Since the solar neutrino mass hierarchy is 
known to be normal, the bipolar conversion in 
the inverted atmospheric hierarchy is essentially 
pair conversion Vf,Vf, — s- Uxi^x from the second 
highest state mi to the lowest state 7713, where 
Vx is a combination of i^^ and Vt- We stress 
that for the inverted hierarchy, bipolar oscilla- 
tions occur for arbitrarily small 6*13. This can 
be used as experimental test by observing a su- 
pernova with mega-ton detectors that are mostly 
sensitive to electron-antineutrinos via the reac- 
tion Ve + p ^ n + [21] : If either the atmo- 
spheric hierarchy is normal or if it is inverted 
with sin^ 6*13 > 10"'^, the bipolar transition is fol- 
lowed by an adiabatic MSW transition, and the 
electron anti- neutrino fluxes seen in two detectors 
behind and in front of the Earth are different. In 
contrast, if the atmospheric hierarchy is inverted 
and sin^ 6*13 < 10^^, the bipolar transition is fol- 
lowed by a non-adiabatic MSW transition, and 



the electron anti-neutrino fluxes seen in two such 
detectors are equal. 

3. Kinematic and Self-Induced Decoher- 
ence between Angular Modes 

The momentum modes p essentially consist of 
energy and angular modes. For spherical symme- 
try one has the radius r as integration variable, 
such that 



P 



- „:„2 , 



{E ~\p\,u = sur Vr 

QE,u,r 1 



(11) 



where Or is the emission angle relative to the ra- 
dial direction at the neutrino sphere, r = R. The 
energy bins will lead to spectral splits |3|7|9|10j . 
whereas angular bins can give rise to kinematic 
and self- induced decoherence |5l6j . To see this, 
we define the flux matrices 

With the radial velocity Vu,r — v-p ^ = cos 6^. = 
■\/l — u{R/r)'^ one has the integral flux and num- 
ber density matrices 



Nr 



r 



d^p 



■ Qp = 



= du 



dEJ, 



du 



dE 



E,u,r ; 



(13) 



/O "'0 '^u,r 

For an average matter velocity v,, one then has 
the general equations of motion: 

VLe + A(r)L(l - V^^rVe) 



(14) 



-V2Gf 



R' 



Nr - Nr 



[Jr Jr) I 1 Je.u.t 



and analogously for antineutrinos with the op- 
posite sign for U,e- Decoherence can occur due 
to the the u— dependence of the velocity Vu,ri ei- 
ther from the flux term in the self-interactions 
(self-induced decoherence) or kinematically in the 
matter terms. The matter term can be trans- 
formed away up to the u— dependence of Vu,r |2I4| . 
The matter-induced multi-angle effect becomes 
important when 



> 



(15) 
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If the matter density is very much larger than 
the neutrino density, the effective oscillation fre- 
quencies of different polarization vectors vary so 
greatly that they stay pinned to the L direction 
and no collective oscillations occur. 

Numerical simulations have been performed to 
determine under which conditions decoherence 
occurs. For ui = 0.3 km~^, sin26' = 10~^ and 
sin26' = 10^'^, the demarcation lines between co- 
herence and decoherence for inverted and normal 
hierarchies, respectively, in the /i-e-plane, can be 
approximated by [6j 



em « 0.225 + 0.027 log^ 
ENH « 0.172 + 0.087 logio 



106 kmT^ ' ' 



106 km" 



-) . (16) 



4. Energy Modes and Spectral Splits 

Integrating Eq. (|14p over E and u and writing 
for the two flavor case = [F^ — Fj^ + T) ■ cr] /2, 
in the limit of small mixing angles one gets [7j 



drB = X M , 



(17) 



where M is an integral of Je.u.t which is not inter- 
esting for our purposes. Thus, ( — Jr)22 — {Jr — 
Jr)ii = F^^ - Fp^ - {F^^ - Fy^\ = const, which 
corresponds to the approximate conservation of 
(P — P) • L in the single mode approximation. 
Since the total lepton number is conserved, elec- 
tron and a;— lepton numbers are conserved sepa- 
rately. 

The spectral split is governed by lepton number 
conservation for both flavors separately: Anti- 
neutrinos swap completely in a bipolar transition. 
To compensate, neutrinos can only swap above a 
certain energy because under typical supernova 
conditions 



> \F. 



(18) 



The spectral splits are created by an adiabatic 
transition between the regime dominated by 
neutrino-self interactions and the low neutrino 
density regime [7] . If this transition is not com- 
pletely adiabatic, as may be the case in a real 
supernova, the spectral split tends to be washed 



out [8]. Spectra splits may be observable in fu- 
ture observations of a galactic supernova explo- 
sion 14J. Such splits may also occur in the an- 
tineutrino sector |10ll3j . albeit at lower energies 
and it is not clear if such features are not washed 
out when taking into account angular modes [13j . 

5. Three-flavor EfTects 




r[km] 



Figure 1. Radial evolution of the fluxes of i^e 
(red), h'fj, (black) and (blue) fluxes, normalized 
to the initial i?e-flux, for a fixed neutrino energy 
E = 20 MeV and an inverted atmospheric hierar- 
chy. For the neutrino parameters we use Amf2 = 
Ar ' 



'sol 



= 7.65 X 10"^ eV% Am 



2.4 X 10-3 eV^ sin^ 



12 = 0.304,1in2 



0.01, 



sin 6123 — 0.4, and a vanishing Dirac phase 6 = 0, 
all consistent with measurements [2D]. The mat- 
ter density proflle A(r) = 4 x 10^ {R/r)^kiai~^ 
was assumed, where the neutrinosphere is at 
R = 10 km. The self-interaction term is taken 
as ^i{r) = 7 X 10^ (^A)V(2 - (i?/r)2)km"\ Af- 
ter bipolar conversion, thick lines represent the 
average fluxes, whereas thin lines signify the en- 
velopes of the fast flux oscillations. 



In normal matter, /x and r leptons appear 
only as virtual states in radiative corrections to 
neutral-current and Vr scattering, causing a 
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r[km] 



Figure 2. Same as Fig.[TJ but for the antineutrino 
fluxes. 



shift AVfj,r — V^GpY^^HB between and i>r, 
where ns is the baryon density. It has the same 
effect on neutrino dispersion as real r leptons with 
an abundance 1221 



Yl 



eff 



3V2G 



(2^)2 
2.7 X 10"^ 



Kf) 




Yn 
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(19) 



where rig was assumed to equal the proton den- 
sity and Yn is the neutron fraction of n^. This 
should be compared to the ordinary MSW poten- 
tial AV = v^Gpi^e^iB where Ye — Ue/nB is the 
electron fraction. 

The potential Eq. (HH) would lead to a MSW 
resonance at density p ~ 3 x lO^gcm"'^, corre- 
sponding to A ~ lO^km"^, provided that the ra- 
dius at which this resonance occurs is beyond the 
radii at which collective oscillations occur. The 
and S'e survival probabilities would then be sensi- 
tive to the 023 angle which governs — Ut mixing 
and such a situation could a priori arise in the ac- 
cretion phase of an iron-core supernova [llj . 

However, the v^ — v-r refractive effect is unlikely 
to play any practical role: At the required high 
matter densities the multi-angle matter effect is 
likely to trigger multi-angle decoherence such that 
the fluxes of the different flavors tend to be max- 



imally equilibrated [12]. 

In Figs. [Hand [2] we show the radial dependence 
of neutrino and antineutrino fluxes for a typical 
example of three-flavor oscillations in single mode 
approximation. The asymmetry parameter e = 
0.25 from Eq. ([9]) is sufficiently large to prevent 
self- induced decoherence according to Eq. ([T6|l . 
Since ^i{r) > A(r)(l - Vu,r) ^ \{r){R/rf /2 out- 
side the synchronized oscillation regime, matter 
induced decoherence should also be negligible ac- 
cording to Eq. (|15p . Only one radial mode was 
thus taken into account in this simulation. After 
the synchronized oscillation phase which lasts un- 
til r ~ 100 km, a bipolar transition occurs which 
lasts until r ~ 300 km. During this bipolar phase, 
the survival probability of electron antineutrinos 
falls off roughly as /i(r)^/^, as discussed in Sect. [21 
whereas the electron neutrino survival probability 
approaches the value e = 0.25, as dictated by ap- 
proximate flavor conservation. If energy modes 
would be included in such a simulation, this 
would result in a spectral split such that the 
flux would not swap with the + Vr fluxes below 
a certain critical energy determined by approxi- 
mate flavor conservation. Furthermore, a. v^ — Vt 
MSW transition would occur at r ~ 100 km, were 
it not for the collective effects that prevent such 
a transition. Finally, an MSW resonance occurs 
at r ~ 2500 km for antineutrinos due to the as- 
sumed inverted hierarchy, after which vacuum os- 
cillations remain. 



6. Open Questions and Conclusions 

We first summarize our main conclusions. Neu- 
trino self-interactions can play a major role in 
the oscillations of neutrinos in core collapse su- 
pernovae. The one-mode approximation is mean- 
while well understood and can be thought of as 
a spinning top. It often leads to a surprisingly 
accurate description of the full problem which in 
the simplest case of spherical symmetry requires 
the introduction of modes both for energy and for 
the direction of a given neutrino trajectory with 
respect to the radial direction. Asymmetry be- 
tween neutrinos and anti-neutrinos leads to spec- 
tral splits in the context of energy modes and, if 
sufficiently large, prevents self-induced decoher- 
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ence of angular modes. Matter effects can also 
lead to decoherence if the charged lepton density 
is larger than the neutrino density. This is also 
likely to mask any significant effects of the second- 
order difference between the i'^ and refractive 
index on resulting neutrino fluxes. If matter den- 
sities are very much larger than neutrino densi- 
ties, the multi-mode flavor polarization vectors 
remain pinned to the direction corresponding to 
flavor eigenstates and no flavor conversion occurs. 

There are still unresolved issues in collective 
neutrino oscillations, including their detailed nu- 
merical description in the absence of spherical 
symmetry, which can be relevant, for example, in 
the presence of hydrodynamic turbulence in the 
background of ordinary matter. It is, for exam- 
ple, currently not completely clear, even concep- 
tually, how to describe damping due to the diflFer- 
ent matter profiles "seen" along different neutrino 
trajectories in this context. Another interesting 
question could be if their could be any significant 
dependence of collective oscillations on the Dirac 
phase b [23] . 
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